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I.  INTRODUCTION 


This  Is  the  FY82  final  report  on  the  program  entitled  "Laser  Remote  Sensing 
of  Atmospheric  Pollutants"  supported  by  the  Air  Force  Engineering  and  Services 
Center  (AFESC).  The  effort  Is  part  of  a  larger  ongoing  program  at  Lincoln 
Laboratory  to  develop  laser  remote  sensing  techniques  for  the  detection  of 
chemical  species  In  the  atmosphere.  Previous  research  for  AFESC  Is  documented 
In  the  FY79,  FY80,  and  FY81  Final  Reports.1-^ 

The  specific  tasks  which  were  conducted  during  FY82  for  this  research 
program  consisted  of  the  following:  (I)  the  Initial  development  of  a 
heterodyne-detection,  differential -absorption  LIDAR  (DIAL)  system,  (2)  the 
development  of  a  computerized  data  acquisition  and  processing  system  for  range- 
resolved,  heterodyne-detection  DIAL  measurements,  (3)  preliminary  range- 
resolved,  heterodyne-detection  DIAL  measurements  of  atmospheric  species,  and 
(4)  a  study  of  the  potential  application  of  the  continuously  tunable  Co:MgF2 
IR  laser  for  remote  sensing  of  airfield  pollutants  and  gaseous  species. 

Each  of  the  tasks  are  described  in  detail  in  the  following  sections. 
Supportive  documentation  Is  Included  in  the  appendices. 

II.  INITIAL  DEVELOPMENT  OF  HETERODYNE-DETECTION  DIAL  SYSTEM 

The  direct -detect Ion  10  um  dual-laser  DIAL  system  as  described  In  Ref.  3 
was  modified  to  permit  operation  with  both  direct  detection  and  heterodyne 
detection.  A  schematic  of  the  DIAL  system  Is  shown  in  Fig.  1  and  a  photograph 
of  the  system  Is  presented  In  Fig.  2.  Two  grating-tuned  nonstablllzed  hybrid- 
TEA  CO2  lasers  provided  pulsed,  tunable,  single-frequency  radiation  near 
10  ym.  The  outputs  from  these  two  lasers  were  directed  out  the  laboratory 
window  and  the  LIDAR  returns  collected  with  a  telescope.  Through  use  of  a  beam 
splitter,  the  Individual  pulsed  LIDAR  returns  could  be  detected  simultaneously 


Photograph  of  DIAL  system. 


both  coherently  (heterodyne  detection)  and  noncoherently  (direct  detection). 
Appropriate  power  normalization,  frequency  monitoring,  and  signal  processing 
were  conducted  as  shown  in  Fig.  1.  This  system  permitted  the  direct  comparison 
of  the  signal -to-noise  ratios  and  the  statistical  characteristics  of  direct  and 
heterodyne  detection  of  the  same  LIDAR  returns. 

LIDAR  returns  were  recorded  from  several  diffuse  targets  and  specular  ret- 
roreflectors  at  ranges  up  to  2.7  km.  A  statistical  analysis  of  these  signals 
was  used  to  quantify  the  relative  merits  of  the  two  detection  techniques.  The 
results  are  given  in  detail  in  Appendices  A  and  B.  A  synopsis  of  the  results 
is  given  below. 

The  heterodyne-detection  DIAL  system  was  found  to  have  a  higher  average 
signal-to-noise  ratio  (SNR)  than  that  of  the  direct-detection  system,  at  the 
expense,  however,  of  less  accuracy  in  the  remote  sensing  measurements  at  closer 
ranges.  These  differences  are  due  to  atmospheric  and  target  speckle  induced 
fluctuations  in  the  LIDAR  returns  and  the  different  ability  of  each  detection 
system  to  compensate  for  these  fluctuations.  The  results  can  be  summarized  as 
fol 1 ows : 

1.  For  the  nonoptimized  DIAL  system  used,  the  average  SNR  of  the 
heterodyne  system  was  a  factor  of  1,000  greater  than  that  of 
the  dl rect -detect ion  system;  these  results  agreed  with  theo¬ 
retical  prediction.  It  is  anticipated,  however,  that  this 
factor  would  be  much  lower  (-10  to  100)  for  an  optimized 
detection  system.* 

2.  The  accuracy  of  the  DIAL  measurements  was  approximately  3%  for 
the  direct -detect ion  system  and  20%  for  the  heterodyne-detection 
system.  These  results  reflect  the  greater  fluctuation  level  of 
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the  LIDAR  returns  measured  with  the  heterodyne  system  due  to 
speckle.  The  di rect -detect ion  system  is  capable  of  speckle  or 
aperture  averaging  of  the  returns,  thereby  reducing  the  magnitude 
of  the  fluctuations. 

3.  The  ability  of  each  detection  technique  to  compensate  for 
the  fluctuations  through  signal  averaging  is  limited  by 
changes  in  the  physical  characteristics  of  the  atmosphere. 

These  changes  were  observed  to  be  the  limiting  factor  in 
determining  the  measurement  accuracy. 

4.  Of  importance  was  the  observation  that  the  accuracy  of  the 
measurements  did  not  improve  as  the  square  root  of  the  number 
of  HOAR  pulse  averaged,  but  saturated  at  a  level  determined 
by  the  observed  atmospheric  drifts.  These  results  agreed  with 
theory,  as  presented  in  Appendix  B. 

5.  The  experimental  results  were  used  along  with  appropriate 
analysis  to  quantitatively  predict  the  accuracy  of  laser  remote 
sensing  for  both  range-resolved  and  path-averaged  measurements. 

An  example  of  such  a  prediction  is  shown  in  Fig.  3,  which 
illustrates  the  trade-offs  of  the  increased  signal  strength 
coupled  with  increased  fluctuations  of  the  heterodyne  system 
compared  to  the  increased  accuracy  but  shorter  detection  ranges 
of  the  di rect -detect ion  system  for  a  given  laser  pulse  energy. 

The  results  give  quantitative  values  for  the  sensitivity  and  accuracy  of 
these  DIAL  measurements  and  provide  Insight  into  the  physical  processes  which 
limit  each  of  the  two  detection  techniques. 
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MINIMUM  RANGE  RESOLVED  CONCENTRATION  (ppb) 


Fig.  3.  Predicted  detection  sensitivity  for  CO2  DIAL  system. 
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III.  DEVELOPMENT  OF  COMPUTERIZED  DATA  ACQUISITION  AND  PROCESSING  SYSTEM 
FOR  HETERODYNE-DETECTION  DIAL  SYSTEM 

A  computerized  data  acquisition  and  processing  system  was  developed 
for  use  with  the  heterodyne-detection  DIAL  system.  Two  systems  were  developed 
for  this  purpose.  The  first  system  was  a  modified  version  of  the  previously 
used  system3  and  consisted  of  a  data  acquisition  system  (computer,  CAMAC 
crate,  and  A/D  converters);  this  system  was  used  to  collect  and  analyze  the 
pulsed  LIDAR  returns.  The  second  system  was  comprised  of  a  computer  inter¬ 
faced  to  a  CAMAC  crate  and  transient  digitizer;  this  system  recorded  the 
transient  waveform  of  the  LIDAR  return s.  Both  of  these  systems  were  capable  of 
analyzing  the  data  in  real  time  at  laser  pulse-repetition-rates  of  up  to  50  to 
100  Hz.  Typical  analysis  included  the  calculation  of  the  ENR  ratio,  the 
statisticalcharacteristics  of  the  heterodyne  and  direct-detection  returns,  and 
the  concentration  of  the  absorbing  species. 

IV.  PRELIMINARY  RANGE-RESOLVED,  HETERODYNE-DETECTION  DIAL 
MEASUREMENTS  OF  ATMOSPHERIC  SPECIES 

Preliminary  range-resolved  heterodyne  LIDAR  measurements  were  performed. 
These  measurements  consisted  of  detecting  the  10.57-pm  P ( 18)  CO2  LIDAR 
returns  backscattered  from  atmospheric  aerosols.  Typical  results  are  shown  in 
Figs.  4  and  5  which  show  the  return  LIDAR  signal  as  a  function  of  time  after 
integrating  128  pulses;  Fig.  4  shows  the  background  signal  when  the  beam  was 
blocked  and  Fig.  5  shows  the  Increased  signal  due  to  the  aerosol  backscattered 
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Range-resolved  LIOAR  signal  as  a  function  of  time  with  laser  beam 
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As  seen  In  Fig.  5,  returns  were  obtained  at  ranges  out  to  approximately 
500  meters  (l.e. ,  3.3  us).  The  returns,  however,  were  observed  to  depend 
strongly  on  the  focus  and  geometrical  overlap  of  the  projected  LIDAR  beam  and 
the  receiver  (telescope)  field  of  view.  Such  an  effect  directly  influences 
the  range  dependence  of  the  LIDAR  returns.  Preliminary  studies  were  made  into 
possible  modification  of  the  transmitter/receiver  optical  design;  by  proper 
alignment  and  focus,  returns  out  to  1200  m  were  obtained,  at  the  expense, 
however,  of  reduced  returns  at  close  in  ranges.  These  studies  pointed  toward 
the  requirement  that  a  common  transmitter/recei ver  telescope  should  be  used  for 
range-resolved  LIDAR  measurements  at  close  ranges  in  order  to  ensure  that  both 
optical  systems  are  focused  at  the  same  range  and  have  maximal  overlap  of  their 
two  fields  of  view.  Further  work  will  be  conducted  in  order  to  maximize  this 
overlap  with  the  present  DIAL  optical  system. 

V.  STUDY  OF  THE  POTENTIAL  APPLICATION  OF  Co:MgF2  LASERS  FOR  REMOTE  SENSING 

The  transition-metal -doped  solid-state  laser,  Co:MgF2,  has  been  developed 
at  Lincoln  Laboratory5  and  is  a  potential  laser  source  for  laser  remote  sen¬ 
sing.  Previous  work  with  the  Co:MgF2  laser  has  resulted  in  the  growth  of 
large  Co:MgF2  crystals,  the  operation  of  a  cw  Co:MqF2  laser,  and  recently, 
high-power  pulsed  operation. 

The  Co:MgF2  laser  is  continuously  tunable  over  a  wide  spectral  region  in 
the  range  of  1.6  to  2.3  urn.  In  addition,  through  the  use  of  frequency  shifting 
techniques,  this  frequency  reqlon  may  be  extended  into  the  visible  as  well  as 
mld-IR.  Figures  6  and  7  show  the  anticipated  tuning  ranges  of  the  Co:MqF2 
laser  superimposed  on  the  transmission  spectrum  of  the  atmosphere.  The 
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Co:MgF£  Is  seen  to  potentially  offer  a  wide  range  of  spectral  coverage  not 
presently  available  with  apy  other  single-laser  source. 

In  order  to  assess  the  potential  usefulness  of  the  CoiMgFg  laser  for 
remote  sensing,  a  pulsed,  Q-switched,  tunable  Co.'MgFg  DIAL  system  has  been 
constructed.  Figure  8  Is  a  photograph  of  this  system.  Preliminary  results 
have  been  obtained  so  far,  and  the  results  are  encouraging. 

At  present,  output  power  In  the  range  of  100  mJ  per  pulse  has  been  obtained 
when  the  laser  Is  free-running,  and  10  to  20  mJ  when  Q-switched.  The  Q-switched 
output  pulse  length  Is  approximately  300  ns  and  the  llnewldth  Is  approximately 
0.15  cm“I.  These  results  were  obtained  with  a  single-element  blrefrlngent 
tuning  filter  and  a  1-mm-thick  uncoated  etalon  In  the  cavity.  However,  frequency 
pulling  effects  were  observed  when  the  laser  was  tuned  away  from  the  peaks  of 
the  Co:MgF2  gain  curve.  The  use  of  a  triple-element  blrefrlnqent  tuning 
filter  and  coated  etalon  should  help  stabilize  the  spectral  output  of  the 
laser. 

Future  work  will  further  Investigate  the  frequency  stability  of  the 
Co  :MgF2  laser  and  test  Its  usefulness  for  the  remote  sensing  of  selected 
species  In  the  atmosphere. 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  presented  In  this  report  have  described  the  capability  of  a 
CO2  DIAL  system  for  the  remote  sensing  of  species  In  the  atmosphere.  The 
differences  between  heterodyne  detection  and  direct  detection  have  been  experi¬ 
mentally  determined,  and  the  accuracy  of  each  system  has  been  shown  experimen¬ 
tally  to  be  dependent  upon  several  factors  Including  atmospheric  effects  and 
target  characteristics.  In  summary,  the  direct-detection  DIAL  system  offers 
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higher  accuracy  and  lower  system  complexity  at  the  expense  of  shorter  detection 
range,  while  the  heterodyne-detection  system  has  greater  detection  range  at  the 
expense  of  lower  accuracy  and  greater  system  complexity. 

A  Co:MgF2  DIAL  system  has  been  developed;  preliminary  results  are  en¬ 
couraging  since  they  indicate  that  this  system  will  prove  useful  for  the  remote 
sensing  and  identification  of  species  In  the  atmosphere.  Further  work  to  assess 
its  capabilities  is  being  conducted. 
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Appendix  A 


The  following  is  a  reprint  of  a  journal  article  published  in  Applied 
Optics,  1  March  1983,  entitled  "Experimental  Comparison  of  Heterodyne  and 
Direct  Detection  for  Pulsed  Differential  Absorption  COg  LIDAR". 
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Experimental  comparison  of  heterodyne  and  direct  detection 
for  pulsed  differential  absorption  C02  lidar 

D.  K.  Killinger,  N.  Menyuk,  and  W.  E.  DeFeo 


A  pulsed  dual-wavelength  dual-C02- laser  differential-absorption  lidar  (DIAL)  system  has  been  developed 
which  permits  simultaneous  heterodyne  and  direct  detection  of  the  same  lidar  returns.  This  system  has 
been  used  to  make  an  experimental  comparison  of  the  SNRs  and  statistical  and  temporal  characteristics  of 
the  DIAL  returns  from  several  topographic  targets.  These  results  were  found  to  be  in  general  agreement 
with  theory  and  were  used  to  quantify  the  relative  merits  of  the  two  detection  techniques.  The  measured 
parameter  values  were  applied  to  an  analytical  treatment  to  predict  system  trade-offs  for  the  remote  sensing 
of  atmospheric  species,  with  application  to  both  path-averaged  and  range-resolved  measurements. 


I.  Introduction 

The  use  of  differential-absorption  lidar  (DIAL)  sys¬ 
tems  with  CO2  laser  sources  has  proven  to  be  a  sensitive 
method  for  achieving  remote  measurements  of  selected 
species  in  the  atmosphere.1  In  general,  most  remote 
sensing  measurements  using  C02  DIAL  have  been  made 
using  either  pulsed  C02  lasers  in  conjunction  with  direct 
(noncoherent)  detection2-4  or  cw5-6  and  Q-switched  cw7 
C02  lasers  with  heterodyne  (coherent)  detection.  The 
advantages  of  using  heterodyne  detection  with  pulsed 
C02  DIAL  have  been  studied  theoretically8-10  but  have 
yet  to  be  verified  experimentally. 

In  this  paper,  an  experimental  comparison  is  made 
of  the  relative  merits  of  heterodyne  and  direct  detection. 
Measurements  were  made  using  a  pulsed  dual-wave- 
length  dual-C02-laser  differential-absorption  lidar 
(DIAL)  system  which  permitted  simultaneous  hetero¬ 
dyne  and  direct  detection  of  the  same  lidar  return.  This 
system  was  utilized  to  obtain  a  direct  experimental 
comparison  of  the  SNRs  and  statistical  and  temporal 
characteristics  of  the  DIAL  returns.  The  results 
quantify  the  increase  in  the  average  detected  lidar 
backscatter  intensity  relative  to  the  average  noise  level 
for  heterodyne  compared  with  direct  detection,  and 
they  provide  a  measure  of  the  relative  accuracy  of  the 
two  detection  techniques  when  applied  to  DIAL  mea¬ 
surements. 
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A  description  of  the  experimental  apparatus  is  given 
in  Sec.  II.  Section  III  presents  the  experimental  results 
which  quantify  the  SNR,  temporal  correlation,  and 
statistical  fluctuation  of  the  lidar  returns;  these  results 
are  shown  to  be  in  reasonable  agreement  with  theory. 
The  parameter  values  determined  in  Sec.  Ill  are  used 
in  an  analytical  treatment  in  Sec.  IV  to  provide  physical 
insight  into  the  relative  trade-offs  between  the  two 
detection  techniques.  Finally,  Sec.  V  presents  an 
overview  of  the  conclusions  and  experimental  results. 

II.  Experimental  Apparatus 

A  schematic  of  the  dual-laser  DIAL  system  is  shown 
in  Fig.  1.  Two  grating  tuned  nons*abilized  hybrid-TEA 
C02  lasers  provide  the  pulsed,  tunable,  single -frequency 
radiation  near  10  n m.  Each  laser  is  similar  to  that  re¬ 
ported  previously1112  but  with  an  additional  36-cm  long 
low-pressure  gain  cell  placed  within  the  1.32-m  cavity; 
the  low-pressure  gain  cell  reduced  the  TEM,*,  output 
linewidth  of  the  laser  from  about  four  longitudinal 
modes  to  a  single  frequency.  The  output  coupler  is  a 
93%  reflectivity  l.ft-m  radius-of-curvature  mirror  placed 
on  a  PZT  mount.  The  TEA  laser  operated  at  a  prf  of 
up  to  100  Hz.  The  output  energy  of  the  hybrid-TEA 
laser  was  ~10  m-l/pulse  with  a  pulse  length  of  200  nsec. 
The  low-pressure  (10  Torr)  gain  cell  was  operated 
continuously  at  slightly  above  threshold,  with  an  output 
power  of  100  mW.  No  active  or  passive  stabilization  of 
the  hybrid-TEA  laser  cavity  was  utilized  except  to 
mount  the  mirror  components  on  large  (15  X  15  X  15 
cm)  aluminum  blocks.  Laser  2  was  fired  a  few  psec 
after  laser  1  to  provide  temporal  separation  of  the  two 
pulses. 

The  frequency  of  each  hybrid  -TEA  laser  was  adjusted 
to  an  offset  frequency  of  ~20  MHz  relative  to  that  of  a 
grating  tuned,  stable.  100-mW  cw  C02  laser  local  os- 
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Pig.  1.  Schematic  of  pulsed  dual-wavelength  dual-CXVlaaer  differential-absorption  tidar  (DIAL)  system  providing  simultaneous  heterodyne 

and  direct  detection  of  the  lidar  returns. 


cillator  (LO);  each  LO  had  a  frequency  stability  of  better 
than  1  MHz.  The  stability  of  the  hybrid -TEA  laser 
compared  to  the  LO  frequency  was  monitored  by  mea¬ 
suring  the  heterodyne  beat  frequency  using  a  room- 
temperature  HgCdTe  detector,13  spectrum  analyzer, 
and  gated  frequency  counter.  A  l-jtsec  gated  frequency 
counter  was  used  to  monitor  the  beat  frequency  of  the 
cw  portion  of  the  hybrid-TEA  laser  radiation  just  prior 
to  the  TEA  pulse;  a  precision  100-nsec  gated  frequency 
counter  was  used  to  measure  the  frequency  of  the  pulsed 
portion  of  the  heterodyne  beat  frequency  during  the 
200-nsec  TEA  pulse.  These  measurements  established 
that  the  frequency  of  the  pulsed  portion  of  the  hybrid- 
TEA  laser  output  tracked  the  frequency  of  the  cw  por¬ 
tion  within  dhl  MHz.  It  was  found  that  the  pulse-to- 
pulse  stability  of  the  hybrid-TEA  laser  frequency  was 
approximately  ±5  MHz  centered  about  the  20-MHz 
offset  frequency,  and  the  intrapulse  frequency  variation 
(chirp)  was  ~1  MHz.  The  TEA  discharge  disrupted  the 
lasing  of  the  cw  gain  cell  because  of  the  formation  of 
contamination  products;  the  time  required  to  reestab¬ 
lish  the  cw  lasing  condition  effectively  limited  the  prf 
to  ~100  Hz,  dependent  on  discharge  and  gas-mixture 
conditions.  It  should  be  added  that,  while  this  sim- 
ple-to-build  hybrid-TEA  laser  does  not  offer  the  sub- 
MHz  stability  required  for  coherent  Doppler  mea¬ 
surements,  its  ±5  MHz  stability  is  more  than  adequate 
for  heterodyne  detection. 

The  outputs  of  the  two  hybrid  lasers  were  joined  at 
a  50/50  beam  splitter.  Portions  of  the  combined  beams 
were  sampled  using  beam  splitters  and  pyroelectric 
detectors  to  measure  the  energy  and  absorption  char¬ 
acteristics  of  each  laser  pulse  against  a  known  gas  in  a 
laboratory  absorption  cell.  The  laser  beam  was  ex¬ 
panded  by  a  X 10  beam  expander  to  ~6-8  cm  in  diam¬ 


eter  and  directed  out  the  laboratory  window  toward 
topographic  targets.  The  backscattered  lidar  radiation 
was  collected  using  an  // 3  Cassegrain  telescope  (variable 
10-30-cm  diam)  and  Bplit  into  two  beams  with  a  50/50 
beam  splitter.  Half  of  the  returned  beam  was  detected 
by  a  cooled  1-mm  diameter  direct-detection  photo- 
conductive  HgCdTe  detector  (NEP  ~4  X  10~9  W),  and 
the  other  half  was  mixed  with  the  C02  local  oscillator 
output  and  detected  by  a  cooled  lOO-^m  diam  1.5-GHz 
heterodyne-detection  photovoltaic  HgCdTe  de¬ 
tector.13 

The  R.F.  heterodyne  signal  was  amplified  using  a 
low-noise  45-dB  gain  amplifier  chain  passed  through  a 
50-MHz  bandwidth  filter,  full-wave  rectified  using  a 
R.F.  transformer  and  zero-bias  square-law  Schottky 
detector  diodes,  envelope  detected,  and  impedance 
matched  with  an  emitter  follower.  The  detector  diodes 
were  operated  in  the  square-law  region,  where  the  out¬ 
put  voltage  is  proportional  to  the  square  of  the  input 
voltage  and  not  in  the  higher  linear  operating  range. 

The  outputs  from  the  HgCdTe  detectors,  frequency 
counters,  and  pyroelectric  detectors  were  sent  to  a  gated 
high-speed  analog-to-digital  data  acquisition  system. 
This  &ystem  monitored  the  beat  frequency  (frequency 
offset)  of  each  laser  pulse  to  ensure  that  it  fell  within  the 
limits  of  the  50-MHz  bandpass  filter,  normalized  the 
individual  lidar  returns  to  the  transmitted  laser  pulse 
energy,  and  calculated  the  statistical  and  temporal 
characteristics  of  the  returns  in  real  time.  In  addition, 
a  transient  digitizer  was  used  in  conjunction  with  an¬ 
other  computer  to  monitor  the  R.F.  lidar  return  wave¬ 
form  for  qualitative  diagnostics.  It  should  be  noted 
that  the  real-time  analysis  of  the  data  acquisition  sys¬ 
tem  (computational  and  graphical  display)  limited  the 
effective  prf  of  the  lidar  system  to  ~10-15  Hz. 
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Fig  1.  CRT  display  showing  dual-laser  lidar  returns  from  a  target 
at  a  range  of  2.7  km  (time  of  flight  of  17  nsec).  Top  trace  is  the  di¬ 
rect-detection  returns,  and  the  bottom  trace  is  the  heterodyne-de- 
tection  returns;  the  transmitted  laser  pulse  at  zero  delay  time  is  slso 
evident  in  the  bottom  trace.  Temporal  separation  between  lasers  1 
and  2  was  5  nsec. 


Figure  2  is  a  photograph  of  a  CRT  display  showing 
the  simultaneous  lidar  returns  using  a  diffuse  target 
(flame-sprayed  aluminum  plate)14  at  a  range  of  2.7  Ion. 
The  upper  trace  shows  the  direct-detection  returns  for 
lidars  1  and  2,  and  the  lower  trace  shows  the  corre¬ 
sponding  envelope-detected  heterodyne  returns.  To 
properly  compare  the  relative  intensity  of  the  returns 
for  the  two  detection  techniques,  it  should  be  noted  that 
the  current  (and  hence  voltage  output)  for  the  direct- 
detection  detector  is  proportional  to  the  square  of  the 
received  optical  electric  field,  Er  Er,  which  is  the  in¬ 
tensity  of  the  received  lidar  return.  For  heterodyne 
detection,  the  output  current  of  the  HgCdTe  detector 
at  the  beat  frequency  is  proportional  to  the  product  of 
E'r  El o,  when  Elo  is  the  optical  electric  field  from  the 
local  oscillator.16  The  square-law  detector  chain  is 
operated  so  that  the  output  voltage  V0  is  proportional 
to  the  square  of  the  input  voltage  Vfa  so  that  Vo  «  Vfo 
<*  E%El o-  Therefore,  the  voltage  signals  seen  in  Fig.  2 
are  proportional  to  the  received  lidar  intensity  for  both 
direct  and  heterodyne  detection. 

m.  Experimental  Data 

A.  Average  SNR  of  Lidar  Returns 

The  average  SNR  of  the  lidar  returns,  which  repre¬ 
sents  the  ratio  of  the  average  intensity  to  the  average 
noise,  was  measured  for  different  targets  at  several 
ranges.  Calibrated  attenuators  were  UMd  in  the  het¬ 


erodyne-detection  optical  path  to  permit  an  accurate 
comparison  of  the  two  detection  techniques.  In  general, 
the  measured  SNR  value  for  the  heterodyne-detected 
returns  was  found  to  be  a  factor  of  ~ 500- 2000  greater 
than  that  for  the  direct-detected  returns.  As  an  ex¬ 
ample,  using  the  side  of  a  large,  painted,  corrugated 
metal  building  at  a  range  of  2.7  km  as  a  diffuse  reflecting 
target,  the  direct-detection  SNR  was  ~10,  and  the 
heterodyne-detection  SNR  after  optical  attenuation  of 
1000  was  ~10  (i.e.,  SNR  ~10,000).  This  observed  dif¬ 
ference  in  the  measured  SNR  values  may  be  compared 
with  theory  (valid  for  SNR  ~1).  The  noise  equivalent 
power  (NEP)  of  the  noncoherent  detector  is  ~4  X  10-9 
W.  The  noise  of  the  heterodyne  detector  with  sufficient 
local  oscillator  power  is  given15  by 

Ps  =*  hvB/ri,  (1) 

where  r\  is  the  quantum  efficiency  of  the  detector,  and 
B  is  the  bandwidth  of  the  detector  amplifier  in  hertz. 
Using  typical  values  for  our  operating  system  of  v  -  2.8 
X  1013  Hz,  rj  ~  0.5,  and  B  =  50  MHz,  Eq.  (1)  indicates 
a  heterodyne-detection  noise  value  of  2  X  10~12  W.  The 
ratio  of  the  estimated  heterodyne-detection  and  the 
direct-detection  noise  value  is  ~2000  and  in  reasonable 
agreement  with  the  observed  experimental  values  for 
the  ratio  of  the  observed  SNR  values  of  1000;  this 
agreement  is  actually  better  than  indicated  since  the 
experimental  SNR  was  ~10.13 

It  is  also  instructive  to  compare  the  observed  lidar 
return  intensity  with  that  predicted  theoretically.  The 
return  lidar  signal  Pr  is  given12  approximately  by 
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DIRECT 
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Fig.  3.  Computer  display  showing  temporal  history  and  statistical 
distribution  (histogram)  of  the  lidar  returns  from  a  diffuse  target 
(metal  building)  at  a  range  of  2.7  km.  The  prf  of  the  lidar  was 
~10  Hi. 
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TaM*  L  Miaaured  Fluctuation  Laval  (Standard  Oavlation)  ot  Lidar 
Returns  From  Targata  at  a  Ranga  at  2.7  km 
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Pr  -  \PtpAK  eip<-2/?R)j'»R2,  (2) 

where  Pr  is  the  transmitted  laser  power,  p  is  the  target 
reflectivity,  A  is  the  telescope  collection  area,  K  is  the 
overall  optical  efficiency,  /S  is  the  absorption  coefficient 
of  the  atmosphere,  and  R  is  the  lidar  range.  Typical 
estimated  values  for  our  lidar  system  are  Pj-  =  l  mJ/100 
nsec  =  104  W,  p  -  0.1,  K  =  0.1,  A  =  0.06  m2,  and  /3  = 
0.125  km-1  at  10.52  pm.  For  a  range  of  2.7  km,  Eq.  (2) 
predicts  a  value  for  Pr  of  1.4  X  10-7  W.  Comparing  this 
value  as  an  estimate  of  the  average  signal  return  with 
the  associated  NEP,  one  estimates  a  SNR  value  of  (1.4 
X  10~7  W)/(4  X  10~9  W)  =  35  for  the  direct-detection 
system  and  (1.4  X  10"7  W)/(2  X  10"12  W)  =  70,000  for 
the  heterodyne-detection  system.  Taking  into  account 
the  order-of-magnitude  estimates  used  in  the  analysis, 
the  calculated  values  for  SNR  are  in  reasonable  agree¬ 
ment  with  the  experimental  values  of  ~10  and  10,000, 
respectively.  It  should  be  noted  that  experiments 
which  investigated  the  effect  of  the  telescope  aperture 
size  on  the  return  intensity  indicated  that  the  hetero¬ 
dyne  signal  increased  by  a  factor  of  ~2  when  the  aper¬ 
ture  was  increased  from  10  to  25  cm  in  diameter  in 
agreement  with  recent  theoretical  prediction.16 

B.  Statistical  Distribution 

The  standard  deviation,  temporal  history,  and  sta¬ 
tistical  distribution  of  the  lidar  returns  were  recorded 
and  differences  observed  and  quantified  between  the 
heterodyne  and  direct-detection  techniques.  Figure 
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Fig.  4.  Temporal  history  and  statistical  distribution  of  the  lidar 
returns  from  a  1-in.  retroreflector  at  a  ranga  of  2.7  km;  prf  »  10  Hz 


3  shows  two  photographs  of  the  graphic  display  of  the 
computerized  data  acquisition  system.  The  upper 
photograph  shows  the  temporal  history  of  the  simul¬ 
taneous  heterodyne  and  direct-detection  10-Hz  prf  lidar 
1  returns  from  the  sides  of  a  large  metal  building  at  a 
range  of  2.7  km.  The  lower  photograph  in  Fig.  3  shows 
the  statistical  distribution  (histogram)  and  related 
parameters  of  the  lidar  returns;  similar  results  were 
observed  for  the  corresponding  return  from  laser  2  and 
also  when  a  1-  X  1-m  flame-sprayed  aluminum  plate  was 
used  as  the  diffuse  target.  As  seen,  the  heterodyne  re¬ 
turns  have  increased  fluctuations  due  to  speckle  com¬ 
pared  with  the  aperture-averaged  direct-detection  re¬ 
turns. 

Figure  4  shows  analogous  results  using  a  1-in.  diam 
retroreflector  as  the  lidar  target  at  a  similar  range. 
Differences  between  these  results  and  those  shown  in 
Fig.  3  are  evident.  It  should  be  noted  that,  for  the  case 
of  the  retroreflector  returns  shown  in  Fig.  4,  attenuators 
were  used  to  reduce  the  intensity  of  the  outgoing  lidar 
beam  of  the  order  of  103-105  so  that  the  relative  am¬ 
plitudes  of  the  signals  shown  in  Figs.  3  and  4  were  ap¬ 
proximately  the  same. 

The  statistical  distribution  of  the  diffuse-target  lidar 
returns  in  Fig.  3  is  seen  to  approximate  the  expected 
negative  exponential  distribution  for  the  heterodyne 
returns  and  appears  to  approximate  a  normal  (Gauss¬ 
ian)  distribution  for  the  direct-detection  returns.17  19 
The  statistical  distribution  of  the  retroreflector  returns 
in  Fig.  4  for  the  heterodyne  case  appears  to  remain  a 
negative  exponential  distribution;  additional  mea¬ 
surements,  however,  have  shown  different  distributions 
under  certain  operating  conditions.20-21  The  distri¬ 
bution  for  the  direct-detection  returns  in  Fig.  4  is  seen 
to  be  similar  to  a  lognormal  distribution. 
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Fig.  5.  Temporal  pulse-pair  cross-correlation  coefficient  as  a  func¬ 
tion  of  pulse  separation  time  between  lasers  1  and  2  for  heterodyne- 
detection  lidar  returns  from  a  diffuse  target  at  a  range  of  2.7  km.  The 
wavelengths  of  the  two  lasers  were  the  same,  the  PVJOi  line  at 
10.522  fjn\. 
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Table  H.  Measured  Cross-Correlation  Coefficient  p  at  Short  Delay  Times 
tor  Dual-Laser  Lldar  Returns  (Target  Range  2.7  km) 


Direc 

t  detec  tion 

Heterodyne  detection 

Target 

x,  =  v 

i  A|  x 

X,  =  X, 

Aj  x  A2 

Diffuse 

0.f» 

o.r> 

0.9 

0.2° 

Retroreflector 

0.9 

0.9 

0.9 

0.9 

a  Dependent  on  target  surface  roughness  and  Ai  —  A2;  the  veiue 
shown  is  for  A[  —  \  >  555  0.016  ^m. 


The  statistical  distributions  shown  in  Figs.  3  and  4  are 
consistent  with  previously  measured  results  for  the  di¬ 
rect-detection  system  as  well  as  those  predicted  from 
theory  and  associated  measurements  for  the  heterodyne 
technique.1-  17  20  The  exact  distributions  shown  in 
Figs.  3  and  4  will  have  an  influence  on  the  accuracy  of 
a  DIAL  measurement  since  they  affect  the  standard 
deviation  of  the  fluctuations.  The  standard  deviation 
will  influence  the  accuracy  of  a  DIAL  measurement 
since  it  determines  the  error  in  the  estimate  of  the  mean 
value  of  the  lidar  returns. 

The  measured  standard  deviation  of  the  lidar  returns 
for  the  two  different  targets  at  a  range  of  2.7  km  is  shown 
in  Table  I.  As  seen,  the  normalized  standard  deviation 
of  the  fluctuation  in  the  lidar  returns  was  found  to  be 
~20%  and  100%  for  the  diffuse  target  returns  and  70% 
and  120%  for  the  1-in.  retroreflector  return  for  direct 
detection  and  heterodyne  detection,  respectively;  re¬ 
sults  obtained  using  the  aluminum  plate  or  the  metal 
building  as  the  diffuse  target  were  essentially  the  same. 
These  values  are  in  order-of-magnitude  agreement  with 
theoretical  predictions-’2  -’ 1  and  quantify  the  differences 
between  the  two  detection  techniques.  Further  analysis 
of  these  results  will  be  presented  in  later  sections. 

C.  Temporal  Correlation 

The  temporal  pulse-pair  cross-correlation  coeffi¬ 
cient24  p  of  the  two  return  pulses  for  the  two  lasers  was 
measured  as  a  function  of  the  temporal  separation  be¬ 
tween  the  two  laser  pulses.  The  correlation  coefficient 
is  a  measure  of  how  well-correlated  the  fluctuations  of 
lidar  1  returns  are  compared  with  the  fluctuations  of 
lidar  2  returns. 

Figure  5  shows  results  obtained  for  the  correlation 
coefficient  of  the  dual-laser  heterodyne  returns  from  a 
diffuse  target  when  the  two  laser  wavelengths  were  the 
same.  As  seen,  p  is  near  unity  at  close  separation  times 
and  falls  off  at  typical  atmospheric  decorrelation  times 
of  the  order  of  10-50  msec.24  The  value  of  p  at  short 
(<50  psec)  delay  times  was  found,  however,  to  be  de¬ 
pendent  on  the  nature  of  the  target,  the  relative  wave¬ 
lengths  of  the  pulse  pairs,  and  the  type  of  detection 
system  used.  Table  II  shows  typical  approximate 
(±0.1)  values  of  p  measured.  It  should  be  added  that, 
while  typical  of  our  experimental  conditions,  the  ob¬ 
served  atmospheric  decorrelation  times  shown  in  Fig. 
5  are  dependent,  in  general,  on  atmospheric  conditions 
(wind  speed,  wind  direction,  and  Cn2  value),  the  di¬ 
mensions  of  the  transmit/receive  optics,  and  the  dy¬ 
namics  of  the  scattered  far-field  speckle  pattern.'’  '24 
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Fig.  6.  Standard  deviation  in  the  estimate  of  the  mean  of  the  lidar 
returns  as  a  function  of  the  number  of  pulses  integrated  for  returns 
from  a  diffuse  target  at  a  range  of  2.7  km. 


The  value  of  p  is  seen  in  Table  II  to  be  large  for  re¬ 
troreflector  returns  for  all  cases,  including  those  in 
which  the  wavelength  of  the  two  lasers  differed  from 
each  other  by  as  much  as  1  pm.  This  result  is  most 
probably  due  to  the  relative  physical  size  of  the  retro- 
reflector  compared  with  the  size  of  the  modulated  in¬ 
tensity  patterns  within  the  transmitted  laser  beam. 
One  estimates  that  under  the  conditions  of  our  experi¬ 
ments,  the  size  of  the  intensity  modulations  (i.e..  co¬ 
herence  length)  was  ~10-50  cm.22  25  In  this  case,  the 
retroreflector  is  essentially  sampling  the  modulated 
intensity  patterns  within  the  transmitted  laser  beam 
and,  as  expected,  should  not  be  a  strong  function  of 
wavelength  or  of  the  detection  technique. 

The  results  for  the  diffuse  target  show,  however, 
strong  dependence  of  the  detection  technique  and  of  the 
laser  wavelength  separation.  While  the  direct-detec¬ 
tion  results  are  consistent  with  previous  results,12-24  the 
heterodyne-detection  results  are  notable  for  two  rea¬ 
sons.  First,  the  high  value  of  p  observed  in  Fig.  5  at 
short  delay  times  indicates  that  the  same  interference 
or  speckle  pattern  is  being  detected  by  both  lidars  1  and 
lidar  2;  the  falloff  of  p  at  time  separations  of  the  order 
of  10-50  msec  indicates  that  the  detected  speckle  pat¬ 
tern  is  being  modulated  by  atmospheric  turbulence. 
Second,  p  is  seen  in  Table  II  to  be  significantly  reduced 
for  heterodyne  detection  from  a  diffuse  target  case  when 
\i  differs  from  X2  by  as  little  as  0.016  pm,  the  difference 
between  two  adjacent  C02  laser  lines.  This  result  is 
most  probably  due  to  the  effect  of  decorrelation  of  the 
speckle  pattern  caused  by  the  surface  roughness  of  the 
target  when  the  wavelength  of  the  radiation  is  varied.21’ 
The  functional  relationship  of  p  on  surface  roughness 
and  wavelength  may  be  expressed  as  o  = 
exp[-(Afc<7j/v2)2],  where  is  the  standard  deviation 
of  the  surface  roughness,  and  Ak  is  the  difference  in  the 
laser  frequencies.26  As  an  order-of-magnitude  esti¬ 
mate,  one  may  take  an  estimated  value  of  o,  =  2  mm 
and  Ak  s  2rrAX/X2  =  0.00091  (pm)-1.  Using  these 
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CORRELATION  TIME  INTERVAL  (s) 

Fig.  7.  Temporal  autocorrelation  coefficient  of  the  lidar  returns  (prf 
*  10  Hz)  as  a  function  of  temporal  separation  between  the 
measurements. 


values,  one  predicts  a  value  of  p  ~  0.2,  in  reason¬ 
able  agreement  with  the  observed  value  shown  in 
Table  II. 

The  above  results  obtained  for  the  value  of  p  for  the 
two  different  detection  techniques  are  of  interest  be¬ 
cause  they  influence  the  degree  to  which  the  effect  of  the 
lidar  return  fluctuations  can  be  reduced  by  the  use  of 
a  dual-laser  system.27  In  particular,  the  low  value  of  p 
observed  for  the  heterodyne  and  diffuse  target  case 
when  \!  j*  X2  with  the  two  wavelengths  separated  by 
as  little  as  0.016  pm  implies  that  little  reduction  will  be 
obtainable  for  such  a  DIAL  system.  The  detailed  im¬ 
plications  of  such  effects  will  not  be  covered  here  but 
will  be  published  in  a  later  paper  which  will  combine  the 
effect  of  both  cross-correlation  and  long-term  atmo¬ 
spheric-induced  effects  on  DIAL  remote-sensing  er¬ 
rors.28 

D.  Measurement  Accuracy  in  Determination  of 
Average  Lidar  Return 

The  accuracy  of  a  laser  remote-sensing  measurement 
is  related  to  the  determination  of  the  average  value  of 
the  lidar  return  signal.  Previous  experiments  using 
direct  detection  have  shown  the  limitation  imposed  on 
the  accuracy  of  such  measurements  by  the  atmo¬ 
spheric-turbulence-induced  temporal  correlation  of  the 
lidar  returns.12  To  establish  this  relative  accuracy  for 
both  heterodyne  and  direct  detection,  the  standard 
deviation  of  the  mean,  tf*,29  of  the  lidar  return  was 
measured.  Figure  6  shows  the  measured  standard  de¬ 
viation  of  the  mean  value  of  the  10-Hz  prf  lidar  returns 
as  a  function  of  the  number  of  pulses  integrated,  n.  As 
seen,  <r„,  and  hence  the  fluctuation-induced  error  of  the 
direct-detection  system,  is  lower  than  that  for  the  het¬ 
erodyne  system,  reflecting  the  decreased  fluctuations 
in  Fig.  3  and  Table  I. 

The  deviation  from  the  expected  square-root  de¬ 
pendence  on  n,  the  number  of  pulses  integrated,  is  due 
to  the  nonzero  temporal  autocorrelation  of  the  returns.30 
Figure  7  shows  the  measured  temporal  autocorrelation 
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Fig.  8.  Calculated  minimum  detectable  path-averaged  concentration 
of  ethylene,  JVmlril  as  a  function  of  range  for  heterodyne-  and  direct- 
detection  DIAL. 


coefficient  for  the  same  data  as  in  Fig.  6  and  clearly  in¬ 
dicates  the  nonrandom  nature  of  the  returns.  It  is  in¬ 
teresting  to  note  that  both  detection  techniques  display 
similar  long-term  correlation  effects  indicative  of  similar 
absorption  or  scattering  properties  of  the  atmosphere 
under  these  conditions.  The  results  shown  in  Fig.  6 
indicate  an  effective  limit30  to  the  improvement  in  the 
SNR  (i.e.,  standard  deviation  of  the  mean)  of  the  order 
of  6  or  7  for  both  the  heterodyne-  and  the  direct-de¬ 
tection  cases  under  the  conditions  of  our  experiments. 
The  implications  of  this  on  laser  remote -sensing  accu¬ 
racy  will  be  discussed  in  Sec.  IV. 

It  should  be  added  that  the  results  shown  in  Figs.  5-7 
are  fairly  typical  results  yet  are  indicative  of  the  par 
ticular  experimental  conditions  (atmosphere,  range, 
system  parameters)  while  the  measurements  were  being 
made.  In  general,  significant  changes  in  the  measured 
values  of  on  and  p  can  occur  over  time  periods  as  short 
as  a  few  seconds  to  a  few  minutes  (see  Figs.  11  and  12  in 
Ref.  29). 

IV.  Comparison  of  DIAL  Measurement  Accuracy  for 
Atmospheric  Species 

The  previous  sections  presented  experimental  results 
which  indicated  the  differences  in  the  signal  strength 
and  fluctuation  levels  of  the  two  detection  techniques. 
In  this  section,  these  results  will  be  used  along  with 
appropriate  analysis  to  numerically  predict  the  accuracy 
of  a  laser  remote-sensing  measurement  for  heterodyne 
detection  and  for  direct  detection. 

Previous  analysis  has  shown  that  the  accuracy  of  a 
DIAL  measurement  may  be  limited  by  the  effect  of  the 
fluctuations  of  the  lidar  returns  at  close-in  ranges  and 
by  the  noise  of  the  detector  at  far  ranges.12  29  To  ex¬ 
press  this  limitation,  one  may  relate  the  minimum 
path -averaged  concentration  of  gas  that  can  be  de¬ 
tected,  Nmjn,  to  the  NEP  of  the  detector  as 
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Fig.  9.  Calculated  minimum  detectable  range-resolved  concentration 
of  ethylene,  Nmin,  aa  a  function  of  range.  Range  resolution  is  500  m. 


Nam  =*  (NEP)xR/{2poaKAPT  e*p(-2j9R)l,  (3) 

where  <r„  is  the  absorption  cross  section  of  the  gas. 

The  corresponding  equation  for  the  fluctuation- 
limited  case  is 

Nain  «*  a„/(2a,R),  (4) 

where  <r„  is  the  standard  deviation  of  the  mean  value  of 
the  lidar  returns. 

Equations  (3)  and  (4)  may  be  used  to  compare  the 
detection  accuracy  of  both  heterodyne-  and  direct- 
detection  DIAL  systems.  As  an  example,  one  may  use 
the  same  parameter  values  as  those  used  earlier  in  Eq. 
(2).  In  addition,  reasonable  values  for  <r„  may  be  ob¬ 
tained  from  Fig.  6;  for  n  =  512,  a„  ■  5%  and  17%  for  the 
direct  and  heterodyne  cases,  respectively.  A  nominal 
value  for  oa  is  32  X  106  (km  atm)-1,  representative  of 
the  absorption  cross  section  of  ethylene  near  10.532 
lim.31 

Using  these  values  in  Eqs.  (3)  and  (4),  one  obtains  the 
curves  shown  in  Fig.  8.  The  envelope  of  detection  for 
each  system  is  the  area  of  the  region  bounded  by  the  two 
detection  limits  for  each  system,  respectively.  As  seen, 
at  close-in  ranges,  the  sensitivity  or  accuracy  of  the  di¬ 
rect-detection  DIAL  system  is  better,  while  at  long 
range,  the  increase  in  average  SNR  of  the  heterodyne 
system  is  responsible  for  the  improved  detectivity  at 
these  ranges.  Figure  8  succinctly  illustrates  the  ap¬ 
propriate  trade-offs  of  the  increased  signal  strength 
coupled  with  increased  fluctuations  of  the  heterodyne 
system  compared  to  the  increased  accuracy  but  shorter 
detection  ranges  of  the  direct-detection  system  for  a 
given  laser  pulse  energy. 

A  similar  analysis  can  also  be  applied  to  range-re¬ 
solved  DIAL  measurements  which  utilize  backscatter 
from  aerosols  in  the  atmosphere  instead  of  topographic 
target  returns.  In  this  case,  for  a  CO2  lidar  system 
which  has  sufficient  sensitivity  to  detect  backscatter 
from  atmospheric  aerosols,  Eqs.  (3)  and  (4)  should  be 


modified  by  replacing  p  with  (/J'A R)  and  oa  with 
(aaAR/R),  where  A R  is  the  range  resolution  of  the 
system,  and  /9'  is  the  Mie  backscatter  coefficient.  For 
the  assumed  values  Pr  ~  100  mJ,  p  ~  0.5  X  10-4  (for  (f 
=  10~7  m-1-sr-1  and  AR  =  500  m),  an  ~  10%  and  2%  for 
heterodyne32  and  direct  detection,33  respectively,  for 
n  -  100,  and  with  the  other  parameter  values  taken  to 
be  the  same  as  in  the  preceding  example,  the  range- 
resolved  detection  curves  shown  in  Fig.  9  are  obtained. 
Figure  9  shows  a  detection  envelope  similar  to  that  in 
Fig.  8,  except  that  the  close-in  detection  sensitivity  is 
independent  of  range  since  Eq.  (4)  has  been  modified 
to  ATmjn  ~  o„/(2<raAR),  and  at  the  long  ranges  Eq.  (3) 
has  been  modified  to  be  proportional  to  (NEP)/?2/ 
exp(-2/ 8R).  In  the  latter  case,  the  relative  increase  in 
detection  range  for  heterodyne  compared  with  direct 
detection  is  dominated  by  the  exp(— 2 0R)  term  at  very 
far  ranges  and  the  1  /R2  contribution  at  closer  ranges; 
this  is  why  the  NEP  curves  for  direct  and  heterodyne 
detection  are  not  related  to  each  other  simply  by  the 
square  root  of  the  ratio  of  the  respective  NEP  values  in 
Fig.  9. 

It  should  be  noted  that  Eqs.  (3)  and  (4)  represent 
approximations  and  were  presented  to  provide  physical 
insight  into  the  relative  detection  limits.  In  the  more 
general  case,  the  values  of  NEP  and  on  as  well  as  the 
range  dependence  in  Eqs.  (3)  and  (4)  would  be  depen¬ 
dent  on  several  factors.  These  include  the  range  de¬ 
pendence  of  the  fluctuation  or  atmospheric  turbulence 
levels,20  aperture-averaging  effects,23  the  use  of  cold 
interference  filters  and  amplifiers  to  reduce  the  di¬ 
rect-detection  NEP  value,13  the  use  of  heterodyne  de¬ 
tector  arrays  to  reduce  speckle-induced  fluctuations,34 
increased  signal  averaging  through  the  use  of  high-prf 
lasers,22  and  the  trade-offs  between  the  effects  of  short- 
and  long-term  temporal  correlation  for  DIAL  remote 
sensing  of  atmospheric  species.28  This  latter  consid¬ 
eration  is  expected  to  be  important  for  range-resolved 
heterodyne  detection  of  atmospheric  aerosols;  as  an 
example,  low  short-term  correlation  values  may  drive 
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the  system  design  toward  the  use  of  high-prf  lasers  in 
order  to  reduce  the  speckle-induced  an  values  under 
certain  operating  conditions.  However,  long-term 
temporal  changes  in  the  atmosphere  may  indicate  the 
advantages  of  a  dual-laser  DIAL  system  in  spite  of  low 
short-term  cross-correlation.28  29 

V.  Conclusions 

This  paper  has  presented  an  overview  of  our  CO2 
heterodyne-  and  direct-detection  DIAL  measurements. 
Using  a  relatively  simple-to-build  hybrid-TEA  CO2 
laser,  lidar  returns  have  been  detected  simultaneously 
both  coherently  and  noncoherently  from  several  targets. 
The  relative  average  SNR  and  fluctuation  levels  of  the 
two  systems  were  measured  and  found  to  qualitatively 
agree  within  factors  of  the  order  of  unity  with  previous 
theoretical  predictions.  Finally,  analytical  equations 
were  presented  which  related  the  detection  sensitivity 
of  the  two  systems  as  a  function  of  range.  The  examples 
presented  clearly  indicate  the  sensitivity  trade-offs  of 
the  two  systems  as  a  function  of  range  and  provide 
physical  insight  into  the  relative  differences.  It  should 
be  noted,  however,  that  while  the  overall  trade-off 
analvsis  as  outlined  in  this  paper  shows  the  relative 
technical  merits,  other  considerations  such  as  system 
complexity  should  be  considered  for  a  complete  system 
performance  analysis. 
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Optics,  15  September  1982,  entitled  "Limitations  of  Signal  Averaging  due  to 
Temporal  Correlation  In  Laser  Remote  Sensing  Measurements". 
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Limitations  of  signal  averaging  due  to  temporal  correlation 
in  laser  remote-sensing  measurements 

N.  Menyuk,  O.  K.  Killinger,  and  C.  R.  Menyuk 


Laser  remote  sensing  involve*  the  measurement  of  laser-beam  transmission  through  the  atmosphere  and 
is  subject  to  uncertainties  caused  by  strong  fluctuations  due  primarily  to  speckle,  glint,  and  atmospheric-tur¬ 
bulence  effects.  These  uncertainties  are  generally  reduced  by  taking  average  values  of  increasing  numbers 
of  measurements.  An  experiment  was  carried  out  to  directly  measure  the  effect  of  signal  averaging  on  back - 
scattered  laser  return  signals  from  a  diffusely  reflecting  target  using  a  direct -detection  differential-absorp 
tion  lidar  (DIAL)  system.  The  improvement  in  accuracy  obtained  by  averaging  over  increasing  numbers 
of  data  points  was  found  to  be  smaller  than  that  predicted  for  independent  measurements.  The  experimen¬ 
tal  results  are  shown  to  be  in  excellent  agreement  with  a  theoretical  analysis  which  considers  the  effect  of 
temporal  correlation.  The  analysis  indicates  that  small  but  long-term  temporal  correlation  severely  limits 
the  improvement  available  through  signal  averaging. 


I.  Introduction 

The  remote  sensing  of  gaseous  species  in  the  atmo¬ 
sphere  using  a  direct-detection  differential-absorption 
lidar  (DIAL)  system  requires  the  measurement  of  the 
average  backscattered  lidar  return  signals  at  both  a 
high-  and  a  low-absorption  frequency  of  the  species 
under  investigation  and  a  determination  of  their  ratio. 
These  measurements  are  subject  to  uncertainties  caused 
by  large  fluctuations  in  the  received  signals  due  pri¬ 
marily  to  speckle,  glint,  and  atmospheric-turbulence 
effects.  The  physics  of  laser-beam  fluctuations  in  the 
atmosphere  has  been  studied  extensively  by  several 
authors,1  with  consideration  given  to  both  one-way  laser 
propagation2-3  and  backscattered  returns  from  remote 
targets.4-8  These  fluctuations  frequently  establish  the 
limit  of  sensitivity  by  being  the  major  cause  of  uncer¬ 
tainty  in  the  average  value  of  the  lidar  signals. 

In  this  paper  we  shall  be  primarily  concerned  with 
determining  the  limitations  that  temporal  correlation 
places  on  the  ability  of  signal  averaging  to  reduce  this 
measurement  uncertainty.  It  should  be  emphasized 
that  this  determination  is  independent  of  the  origin  of 
the  fluctuations  causing  the  uncertainty. 
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The  measurement  uncertainty  is  defined  as  the  nor¬ 
malized  standard  deviation  of  the  lidar  signals.  A 
standard  technique  for  improving  measurement  accu¬ 
racy  is  to  increase  the  number  of  measurements  and 
take  the  average  value.  For  N  measurements,  the  un¬ 
certainty  is  expected  to  decrease  as  N~in,  assuming  that 
all  the  measurements  are  independent.  An  experiment 
was  conducted  which  directly  tested  the  validity  of  this 
assumption  by  measuring  the  standard  deviation  of 
lidar  returns  as  a  function  of  the  number  of  pulses  av¬ 
eraged.  The  reduction  of  the  standard  deviation  by 
signal  averaging  was  found  to  be  much  more  weakly 
dependent  on  N  than  N~1/2,  indicating  that  the  indi¬ 
vidual  lidar  measurements  are  not  mutually  indepen¬ 
dent.  This  lack  of  independence  was  established  by  an 
analysis  of  the  data,  which  showed  the  presence  of 
positive  temporal  correlation  extending  over  time  pe¬ 
riods  of  several  seconds  to  minutes  for  different  data 
samples.  This  long-term  correlation  of  the  lidar  returns 
is  presumably  due  to  slow  changes  in  atmospheric  ab¬ 
sorption  occurring  during  the  measurement  interval. 
Among  the  potential  sources  of  such  changes,  including 
pressure,  temperature,  and  aerosol  variations,  humidity 
fluctuations  are  believed  to  play  the  most  significant 
role. 

The  experimentally  observed  variation  of  the  stan¬ 
dard  deviation  with  the  number  of  pulses  averaged  is 
shown  to  be  in  excellent  agreement  with  a  theoretical 
analysis  which  considers  the  temporal  correlation  of 
successive  lidar  returns.  It  should  be  noted  that  this 
theoretical  analysis  quantifies  the  deviation  of  any 
signal-averaging  process  from  the  N~U2  behavior  due 
to  nonrandom  behavior  and  is,  in  general,  applicable  to 
a  wide  range  of  signal-averaging  phenomena. 
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N.  hpotwwtol  Results 

The  direct-detection  dual-laser  DIAL  system  used 
for  this  experiment  has  been  described  previously.7'8  It 
employs  two  line-tunable  mini-TEA  COj  lasers*  to 
provide  the  pulsed  coherent  radiation.  The  time  delay 
between  the  firing  of  the  two  lasers  was  maintained  at 
50  Msec,  so  that  the  atmosphere  may  be  considered 
frozen  between  firings.7  A  portion  of  each  beam  was 
sent  to  a  pyroelectric  detector  to  normalize  the  energy 
in  each  pulse,  but  the  bulk  of  the  radiation  was  directed 
outside  the  laboratory  toward  a  1-  X  1-m  flame-sprayed 
aluminum  plate  located  at  a  range  of  2.7  km  which 
served  as  a  diffusely  reflecting  target  The  lidar  returns 
were  recorded  and  normalized  in  a  computerized 
data-acquisition  system. 

For  this  experiment  laser  1  was  fired  on  the  10.67-jim 
P(28)  CO2  laser  transition,  and  50  ftaec  later  laser  2  was 
fired  on  the  10.61  -jim  P(22)  transition.  At  these  tran¬ 
sition  frequencies  the  atmospheric-transmission  levels 
are  high,  with  water  vapor  the  primary  source  of  atmo¬ 
spheric  absorption.10 

A  total  of  22,528  normalized  lidar  return  pulses  from 
each  laser  was  recorded  for  later  statistical  sinalysis;  the 
lidar  return  signals  were  normalized  on  a  pulse-to- pulse 
basis  to  the  energy  in  each  transmitted  pulse.  The 
entire  process  took  40  min,  corresponding  to  an  overall 
pulse  repetition  frequency  of  slightly  under  10  Hz. 
Computational  constraints  limited  analysis  to  sets  of 
12,288  normalized  pulse  returns  from  each  laser  and 
their  ratios.  The  initial  and  final  sets  of  12,288  pulses 
each  exhibited  somewhat  differing  behavior.  There¬ 
fore,  a  full  analysis  was  carried  out  for  these  two  sets  of 
data  separately.  It  should  be  noted  that  there  is  an 
overlap  in  the  data  of  the  two  sets;  that  is,  the  last  2048 
pulses  of  the  initial  set  correspond  to  the  first  2048 
pulses  of  the  final  set.  With  T  defined  as  the  total 
number  of  pulses  in  a  data  set  (in  the  present  case  V  » 
12,288),  statistical  analysis  of  the  two  sets  of  lidar  return 
data  included  a  determination  of  the  average  value  of 
each  set  of  T  returns,  the  standard  deviation  of  the  re¬ 
turns,  and  the  normalized  standard  deviation  of  the 


average  value  of  the  returns,  as  as  a  function  of  N,  the 
number  of  pulses  being  averaged. 

The  standard  deviation  was  determined  using  two 
distinct  methods  to  define  the  mean  values  of  the  N 
pulses  being  integrated.  The  first  method  used  a  seg¬ 
mental-averaging  approach.  In  this  case,  as  an  exam¬ 
ple,  for  N  *  8,  the  segments  averaged  would  be  1-8, 
9-16,  17-24,  ....  12,281-12,288,  and  the  measured 
standard  deviation  of  the  mean  is  based  on  the  resultant 
12,288/A/  average  values.  A  second  method  used  a 
running  average  over  the  set,  i.e.,  when  N«8,  the 
groups  averaged  are  1-8,  2-9, 3-10, . . . ,  12,281-12,288, 
and  the  measured  standard  deviation  is  then  based  on 
the  12, 288- A/  average  values.  It  will  be  established  in 
the  course  of  this  paper  that  the  two  methods  yield  es¬ 
sentially  the  same  results  for  a  sufficiently  large  data 
sample. 

A  comparison  of  the  results  obtained  by  the  two 
methods  for  the  initial  set  of  12,288  pulses  from  both 
laser  1,  Lit  and  laser  2,  L2,  is  given  in  Table  I,  which 
shows  the  measured  percentage  signal-averaged  stan¬ 
dard  deviation  <tn  of  the  normalized  lidar  returns  as  N 
is  increased  from  1  to  2048  by  factors  of  2.  Little  dif¬ 
ferences  between  the  two  averaging  methods  is  observed 
for  all  values  of  N,  with  the  results  essentially  identical 
for  N  <  256.  The  corresponding  results  for  the  final  set 
of  12,288  pulse  returns  are  given  in  Table  II.  The  tables 
also  contain  the  results  for  the  ratio  of  the  pulse  pairs; 
the  ratio  values  are  required  for  differential-absorption 
lidar  measurements  (see  Appendix  A).  There  is  no 
observable  reduction  in  the  standard  deviation  using 
the  ratios  in  the  initial  set,  which  is  consistent  with  the 
fact  that  the  measured  cross-correlation  coefficient  of 
the  corresponding  pulses  from  the  two  lasers  was 
<0.5.11,12  However,  the  use  of  ratios  in  the  final  set  is 
seen  to  lead  to  a  marked  reduction  in  the  standard  de¬ 
viation  of  the  mean  at  large  N.  This  improvement, 
which  occurs  despite  a  pulse-pair  cross-correlation 
coefficient  <0.5  (as  seen  by  the  lack  of  improvement  at 
N  «  1),  has  been  shown  to  be  due  to  a  slow  drift  in  at¬ 
mospheric  absorption  which  occurred  during  the  time 
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N 

Laser  1  P(28) 

Laser  2  P(22) 

Ratio  L1/L2 

Running 

average 

Segment*) 

average 

Calculated 

Running 

average 

Segmental 

average 

Calculated 

Running 

average 

Segmental 

average 

Calculated 

1 

20.5 

20.5 

17.7 

17.7 

22.3 

22.3 

2 

16.7 

16.7 

16.7 

14.7 

14.6 

14.7 

17.0 

17.0 

17.0 

4 

13.7 

13.8 

13.7 

12.3 

12.2 

12.3 

13.3 

13.3 

13.2 

8 

11.2 

11.3 

11.2 

10.2 

10.1 

10.3 

10.6 

10.8 

10.7 

16 

9.0 

9.0 

9.1 

8.4 

8.3 

8.6 

8.8 

8.7 

8.9 

32 

7.2 

7.3 

7.3 

7.0 

7.0 

7.2 

7.5 

7.6 

7.7 

64 

5.8 

6.0 

5.8 

6.0 

6.1 

6.2 

6.6 

6.6 

6.8 

128 

4.6 

4.8 

4.8 

5.2 

6.2 

5.4 

5.8 

5.8 

6.2 

256 

3.7 

3.8 

4.0 

4.5 

4.5 

4.7 

4.9 

4.9 

5.2 

512 

3.0 

3.0 

3.0 

3.7 

3.9 

4.0 

3.8 

4.0 

4.3 

1024 

2.6 

2.7 

2.3 

3.0 

3.3 

3.2 

2.8 

3.3 

3.3 

2048 

2.3 

2.2 

2.6 

2.3 

2.8 

2.4 

2.2 

2.4 

2.6 
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N 

Later  1  P(  28) 

Laaer  2  P'-2) 

Ratio  LI/L2 

Running 

average 

Segmental 

average 

Calculated 

Running 

average 

Segmental 

average 

Calculated 

Running 

average 

Segment*) 

average 

Calculated 

1 

21.7 

21.7 

18.3 

18.3 

22.2 

22.2 

2 

18.0 

18.0 

18.0 

15.2 

15.2 

15.2 

16.8 

16.8 

16.8 

4 

15.1 

15.1 

15.1 

12.7 

12.7 

12.8 

12.9 

12.9 

12.9 

8 

12.8 

12.8 

12.9 

10.7 

10.7 

10.9 

10.0 

10.0 

10.0 

16 

10.9 

10.8 

11.1 

9.1 

9.1 

9.3 

7.9 

7.9 

8.1 

32 

9.5 

9.5 

9.7 

7.9 

7.8 

8.0 

6.6 

6.4 

6.8 

64 

8.5 

8.5 

8.6 

7.1 

7.1 

7.1 

5.4 

5.3 

5.7 

128 

7.9 

7.9 

8.0 

6.6 

6.6 

6.6 

4.5 

4.5 

4.8 

256 

7.5 

7.5 

7.6 

6.2 

6.0 

6.2 

3.9 

3.9 

4.1 

512 

7.2 

7.3 

7.2 

5.8 

5.8 

5.9 

3.3 

3.4 

3.5 

1024 

6.9 

7.0 

6.9 

5.5 

5.2 

5.6 

2.7 

2.8 

3.1 

2048 

6.3 

6.8 

6.5 

4.9 

4.8 

5.3 

2.0 

2.4 

2.3 

Fig.  1.  Measured  percentage  standard  deviation  of  lidar  return 
signals  as  a  function  of  the  number  of  pulses  averaged:  (a)  Returns 
from  laser  1  and  laser  2  for  the  initial  set  of  12,288  pulses;  (b)  returns 
from  laser  1,  laser  2,  and  their  ratio  for  the  final  set  of  12,288  pulses. 
The  ratio  results  for  the  initial  set  are  not  shown  in  (a)  as  they  are 
essentially  the  same  as  those  obtained  for  the  individual  lasers. 


interval  encompassed  by  the  final  set.12  The  relatively 
poor  pulse-pair  cross-correlation  obtained  from  the  two 
lasers,  despite  the  effectively  frozen  atmosphere  during 
the  50-peec  interval  between  pulses,  is  attributable  to 
instrumental  effects,  specifically  a  combination  of 
electronic  effects,  pulse-to-pulse  changes  in  mode 
quality,  and  small  directional  shifts  in  the  laser 
beams.7 

The  variation  of  the  segmental  standard  deviation  <7n 
as  a  function  of  N  is  shown  graphically  in  Fig.  1  for  the 
initial  and  final  sets,  respectively.  It  is  seen  that  in  all 
cases  the  reduction  of  as  with  N  is  significantly  smaller 
than  the  N~in  dependence  predicted  for  independent 
measurements. 


III.  Theory 

The  departure  from  N~i/2  behavior  seen  in  Fig.  1 
indicates  that  the  measurements  are  not  truly  inde¬ 
pendent;  that  is,  sucessive  pulses  are  correlated.  This 
behavior  can  be  explained  quantitatively  by  a  theoret¬ 
ical  analysis  which  considers  the  effect  of  temporal 
correlation. 

To  evaluate  the  temporal  correlation  of  the  lidar  pulse 
returns,  let  us  define  /*  =  /(£*)  as  the  normalized  de¬ 
viation  of  the  feth  pulse  return  (occurring  at  time  4) 
from  its  mean  value  P  over  the  full  set  of  data.1 
Then, 


where  P*  is  the  feth  lidar  pulse  return  signal.  The 
normalized  variance  of  the  full  set  of  individual  pulses 
is  defined  as 

a2 = «m2> - ^  i :  H’  (2) 

r  *-i 

where  T  is  the  total  number  of  pulses  in  the  set,  which 
is  equal  to  12,288  in  our  experiments.  The  temporal 
autocorrelation  coefficient  Pj  for  a  delay  time  equal  to 
jr  is  then  given  by 

Pi  -  ~  </«*)/<«*  +  jr))  -  ■  ~~  z  /*/*♦;,  (3) 

where  r  is  the  time  interval  between  pulses  These 
definitions  of  pj  and  a2  are  independent  of  the  proba¬ 
bility-distribution  function  of  the  signals. 

Equation  (3)  was  used  to  calculate  the  correlation 
coefficient  pj  for  j  *  1,2,4, 8, . . . ,  1024,2048,  using  the 
lidar  pulse  returns  from  L  t,  L2,  and  L\tLi-  The  results 
are  given  in  Fig.  2  for  the  initial  and  final  sets  of  data, 
respectively.  It  is  seen  from  the  figures  that  small 
positive  temporal  correlations  persist  in  the  first  set  out 
of  10  sec  and  beyond,  while  the  presence  of  a  slow  drift 
in  the  atmospheric  absorption  during  the  final  set12  is 
seen  to  result  in  a  larger  positive  correlation  which 
persists  beyond  200  sec  in  the  lidar  returns. 
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Previous  studies  of  the  temporal  correlation  coeffi¬ 
cient  for  backscattered  CO2  laser  radiation  through  the 
atmosphere  from  a  hard  target  have  indicated  that  the 
atmosphere  is  effectively  frozen  for  the  order  of  1-5 
msec,  with  a  significant  drop  in  correlation  by  50-100 
msec.7  However,  residual  positive  correlation  has 
been  observed  to  persist  out  to  several  seconds.6 

The  short-term  correlation  is  primarily  due  to  tur¬ 
bulence  caused  by  thermal  fluctuations  in  the  atmo¬ 
sphere.14  For  temporal  correlations  over  the  longer 
time  intervals  considered  in  this  paper,  changes  in  ab¬ 
sorption  due  to  humidity  fluctuations  appear  to  play  a 
significant  role.1-  The  resultant  temporal  autocorre¬ 
lation  as  given  in  Fig.  2  represents  the  summed  effect 
of  all  the  time  constants  involved,  including  those 
mentioned  above. 

To  establish  if  the  temporal  correlation  levels  shown 
in  Fig.  2  are  sufficient  to  cause  the  strong  reduction  in 
the  effectiveness  of  pulse  averaging  relative  to  N~l/ 2 
behavior  as  shown  in  Fig.  1,  the  relationship  between 
<t/v  and  Pj  was  derived  for  both  the  segmental-  and 
running-average  approaches.  The  theoretical  analysis 
for  segmental  averaging  is  given  in  Appendix  B  and 
leads  to  the  relationship 

a  [  v-i  I1/2 

<XN(seg>  =  —  |l  +  2  E(1  -y/Dd  -j/N)pj |  •  (4) 

Using  a  running  average  of  the  lidar  returns,  the  re¬ 
lationship  is  shown  in  Appendix  C  to  be 

a  I  N-\  li/2 

ff*(run)  =  —  1  +  2  E  U-y/A0p,  •  (5) 

VM  /-l  J 

It  should  be  noted  that  in  the  absence  of  correlation 
(p>  =  0),  Eqs.  (4)  and  (5)  predict  the  expected  N~l/2 


TIME  DELAY  (•) 


Fig.  2.  Autocorrelation  coefficient  of  lidar  return  signals  from  laser 
1,  laser  2,  and  their  ratio  as  a  function  of  time  delay  between  signals 
based  on  a  10-Hz  pulse-repetition  rate:  (a)  Initial  set  of  12,288  pulses; 
(b)  final  set  of  12,288  pulses. 


dependence.  It  is  shown  in  Appendix  B  that  Eq.  (5)  is 
an  excellent  approximation  to  Eq.  (4),  so  this  simpler 
form  can  be  used  for  both  segmental-  and  running- 
average  calculations.  This  effective  near-equivalence 
of  a\(seg)  and  tr^(run)  accords  with  the  agreement 
observed  experimentally,  as  seen  in  Tables  I  and  II. 

Values  for  <r,v  were  calculated  on  the  basis  of  Eq.  (5), 
using  the  values  of  p;  given  in  Fig.  2  for  j  =  1,2, 4, 8, .... 
1024,2048,  and  assuming  a  linear  interpolation  for  all 
other  values  of  j.  The  results  of  the  calculation  are 
included  in  Tables  1  and  II.  The  agreement  between 
the  calculated  and  the  measured  values  of  c  s  is  seen  to 
be  excellent  even  at  large  N  values,  where  the  assump¬ 
tions  of  equal  weighting  (discussed  in  Appendices  B  and 
C)  and  the  linear  interpolation  of  p;  values  are  weakest. 
This  agreement  effect  ively  validates  these  assumptions 
and  shows  that  the  presence  of  small  but  persistent 
temporal  correlation  severely  limits  the  improvement 
available  through  signal  averaging  of  lidar  returns  from 
hard  targets. 

It  should  be  noted  that  the  process  considered  in  the 
above  analysis  is  almost  certainly  not  stationary  on  the 
time  scale  of  the  experiments  and  is  definitely  not 
Markovian,  since  the  autocorrelation  function  always 
decays  exponentially  in  Markovian  processes.15-16  Most 
theoretical  work  is  concerned  strictly  with  stationary 
processes,17  and  much  of  it  is  concerned  with  Markovian 
processes.  The  work  presented  here  is  subject  to  nei¬ 
ther  of  these  limitations. 


IV.  Data  Analysis 

To  establish  the  relative  contribution  to  trN  of  the 
short-term  and  the  long-term  autocorrelation  coeffi¬ 
cients,  Eq.  (5)  has  been  evaluated  assuming  (1)  p;  equals 
a  constant,  p;  =  p;  (2)  p,  decreases  logarithmically  with 
increasing  j,  pj  =  A  -  B  ln(j);  and  (3)  pj  decreases  ex¬ 
ponentially,  Pj  =  exp(—  jot).  The  evaluation  for  each 
of  these  three  assumed  functional  relationships  of  p;  is 
given  in  Appendix  D. 

For  the  case  of  a  constant  value  pj  =  p,  Eq.  (D2)  shows 
that  in  the  limit  for  large  values  of  N 

<TN  «  <7\/p  (61 

This  situation  is  approximated  experimentally  when 
there  is  a  continuous  atmospheric  drift  throughout  the 
measurement  period,  as  was  the  case  during  the  final  set 
of  data  taken  with  lasers  1  and  2  as  seen  in  Fig.  1(b).  As 
an  example,  using  a  =  0.22  from  Fig.  1(b)  and  p  =  0.05 
from  Fig.  2(b),  substitution  into  Eq.  (6)  leads  to  as  ~ 
0.05.  This  result  is  consistent  with  the  experimental 
data  shown  in  Fig.  1(b)  and  indicates  that  improvement 
in  the  standard  deviation  of  the  mean  is  limited  in  this 
case  to  a  factor  of  ~4,  independent  of  the  number  of 
pulses  averaged. 

A  more  accurate  approximation  to  the  experimental 
results  can  be  obtained  by  noting  in  Fig.  2  that  there  is 
an  initial  rapid  decrease  in  p;  followed  by  a  p,  value 
which  decreases  slowly  with  increasing  j.  This  behavior 
can  be  approximated  reasonably  by  a  model  which  as¬ 
sumes  pj  to  be  the  sum  of  an  exponential  term  and  a 


3380  APPLIED  OPTICS  /  Vol.  21,  No.  18  /  15  September  1982 

30 


logarithmically  decreasing  term  as  described  in  Ap¬ 
pendix  D. 

A  value  of  a  is  chosen  such  that  the  exponential  term 
p j  =  exp(—  ja)  is  the  dominant  contributor  to  as  at 
small  values  of  j,  but  its  contribution  becomes  negligible 
by  j  >  16.  Beyond  this  value,  the  effect  of  the  expo¬ 
nential  term  alone  on  as  would  be  approximated  by  Eq. 
(D7),  where  it  is  shown  that  as  at  the  larger  values  of 
N  will  be  increased  by  the  constant  factor  (coth(a/2)j1/2 
over  the  noncorrelated  case  but  exhibits  an  N~1/2  de¬ 
pendence.  The  failure  of  as  to  achieve  this  N~iJ7  de¬ 
pendence  is  due  to  the  presence  of  the  logarithmically 
decreasing  term  pj  =  A  —  B  ln(/).  For  the  larger  values 
of  N  in  Fig.  1,  the  resulting  as  values  are  then  almost 
entirely  due  to  a  term  of  this  form. 

These  effects  are  illustrated  in  Fig.  3,  which  shows  the 
values  of  as  obtained  by  applying  Eq.  (5)  for  pj  = 
exp  (-ja),  pj  =  A  -  B  ln(j),  and  p,  *  exp(-;a)  +  A  - 
B  lnO),  where  a  -  0.693,  A  =  0.14,  and  B  =  0.02.  The 
values  of  A  and  B  were  chosen  to  reasonably  approxi¬ 
mate  the  values  of  pj  at  high  j  shown  for  laser  2  in  Fig. 
2(a).  The  value  of  a  was  chosen  such  that  exp(-a)  = 
0.5,  which  approximates  the  experimental  data.  The 
resulting  as  curves  in  Fig.  3  are  of  the  proper  magnitude 
compared  with  the  experimental  results  in  Fig.  2  and 
show  that  the  standard  deviation  after  signal  averaging 
over  a  large  number  of  pulses  is  primarily  due  to  the 
long-term  correlation  from  A  —  B  lnO). 

In  determining  the  extent  to  which  remote-sensing 
accuracy  can  be  improved  by  increasing  the  number  of 
data  points,  it  should  be  noted  that  if  p  is  any  mono- 
tonically  decreasing  or  constant  function  of  delay 
time. 


on  ^  i t 


(7) 


This  is  shown  in  Appendix  D  to  be  true  for  a  constant 
correlation  coefficient  (i.e.,  pj  =  p  —  ps-i)-  Therefore, 
to  prove  the  statement,  it  is  only  necessary  to  note 
that 


a  N- 1  1/2 

"N  *  —r~  1  +  2  I  (1  -  j/N)pj 

viV  ,.i 


O  1 

7*,  +  2£(1 


)!N)ps-\ 


(8) 


since,  by  definition  of  the  correlation  coefficient  as  a 
monotonically  decreasing  or  constant  function  of  delay 
time,  p ,  >  ps-t  for  all  j  ^  N  —  l. 

Equation  (7)  represents  a  fundamental  inequality 
which  can  impose  a  severe  limitation  on  the  improve¬ 
ment  attainable  by  averaging  over  a  large  sample.  As 
long  as  the  monotonically  decreasing  temporal  corre¬ 
lation  coefficient,  ps-\,  has  a  finite  value  for  a  given 
time  interval,  that  value  limits  the  improvement  in  the 
standard  deviation  that  can  be  obtained  by  signal  av¬ 
eraging  regardless  of  the  number  of  pulses  averaged 
during  that  interval.  For  example,  in  the  time  interval 
required  to  bring  the  temporal  correlation  coefficient 
down  to  1%,  it  is  impossible  for  signal  averaging  to 
achieve  more  than  a  tenfold  reduction  in  the  standard 
deviation.  It  should  further  be  noted  that  this  degree 
of  improvement  occurs  only  for  a  constant  correlation 


coefficient.  For  a  decreasing  p, ,  the  improvement  can 
be  significantly  smaller. 

It  is  apparent  from  the  above  considerations  that  very 
small  temporal  correlation  values  can  severely  limit  the 
improvement  achievable  by  signal  averaging.  This  can 
be  an  important  factor  since  our  experimental  data  in¬ 
dicates  that  such  small  correlation  values  may  persist 
over  long  time  intervals,  at  least  for  lidar  returns  from 
a  stationary  target. 

The  Lincoln  Laboratory  portion  of  this  work  was 
supported  by  the  National  Aeronautical  and  Space 
Administration  and  the  Air  Force  Engineering  and 
Services  Center.  The  University  of  Maryland  portion 
of  this  work  was  supported  by  the  Department  of  En 
ergy. 


Appendix  A.  Application  to  DIAL  Measurements 

The  use  of  DIAL  measurements  to  determine  the 
concentration  of  gaseous  species  in  the  atmosphere  re¬ 
quires  determination  of  the  ratio  of  the  average  value 
of  lidar  returns  P  and  P'  at  two  different  wavelengths. 
The  accuracy  in  the  measurement  of  this  ratio,  P/P',  has 
been  shown11  to  be  related  approximately  to  the  accu¬ 
racy  of  each  single  lidar  measurement  as 

<d>//v  =  o%  +  op  —  2poopOp,  (All 

where  <rp/p-,  oP,  and  <r/>-  are  the  normalized  standard 
deviations  of  the  measurement  of  P/P’,  P,  and  P' ,  re¬ 
spectively,  and  pn  is  the  temporal  pulse-pair  cross-cor¬ 
relation  coefficient  of  the  individual  P*  and  Pk  lidar 
measurements.8  "  As  seen  in  Eq.  (A1 ),  if  p0  =  0,  then 
ap/p-  =  aj>+  o'p-;  if  p0  ~  1,  then  aP/P  is  significantly  re 
duced  in  value.  As  ap/p-  is  a  measure  of  tho  accuracy 
of  a  DIAL  measurement,  the  presence  of  a  large  tern 
poral  cross  correlation  leads  to  a  significant  improve¬ 
ment  in  this  measurement  accuracy  relative  to  DIAI . 
experiments  obtained  with  a  system  in  which  the  mea¬ 
surement  at  the  two  frequencies  is  not  made  within  the 
correlation  time. 

Since  the  effect  of  signal  averaging  is  to  operate  on 
both  sides  of  Eq.  (Al)  with  (1  /T)f  dt,  these  results  may 
also  be  applied  directly  to  signal-averaged  DIAL  mea 
surements.  In  this  case,  aP/P',  aP,  and  <xp-  in  Eq.  (Al) 
become  the  values  of  the  signal-averaged  standard  de¬ 
viation  of  the  mean  as  calculated  in  Eqs.  (4)  and  (5),  an>' 
p o  is  the  temporal  cross  correlation  coefficient  of  the 
signal-averaged  values.  An  example  of  the  application 
of  Eq.  (Al)  using  signal-averaged  values  was  presented 
in  Ref.  12,  where  excellent  agreement  between  calcu¬ 
lated  results  and  experimental  data  was  obtained. 


Appendix  B.  Segmental  Averaging 

The  variance  of  the  T  lidar  returns  which  have  been 
segmentally  averaged  over  N  returns  is  defined  by 

i  1  f//i  +  h  +  •  ■  +  In\2  .  //n+i  ■+  •  +  lw\2 

*  m  11  N  )  + 1  N - ) 


Ih-tN-u  + 

■  +  /ry 

l  N 

1 

Expanding  Eq.  (Bl)  leads  directly  to 
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IN  -  2\ 

+  (—|(I-2^+  - 
+  I|r-(N-i)pJv.1)p. 


^  =  —  |l  +  2  E  (1  -;7N)(1  -  j/Dpj 


(B3) 


(B4) 


M  (Number  of  PuIm«  Ay*rog#d) 

Fig.  3.  Calculated  percentage  standard  deviation  as  a  function  of 
the  number  of  pulses  averaged  for  p,  -  exp(-;a),  p,  «  A  -  B  ln(j). 
and  their  sum,  where  a  =  0.693,  A  =  0.14,  B  =  0.02,  and  a  *  0.2. 


The  assumption  of  equal  average  weighting  for  the 
interpolated  terms  appears  reasonable  given  the  well- 
behaved  temporal  correlation  behavior  shown  in  Fig. 
2.  The  ultimate  validation  for  this  assumption  lies  in 


a%  -  TTT;  |(f{  +  /?  +  ...  +  /?■>  +  2  |(/ 1/2  +  /2/;i  +  .  •  +  In-  In 1 
Nl 

+  Un+iIn+2  +  +  IlN-lIlN)  +  ...  +  (/r-(N-l>/|'-(tV-2l  +  +  / 1'  -  l/  rl 

+  2t(/|/3  +  I  Jt  +  . . .  +  /jv-2/n!  +  (//V+1//V+3  +  ■  +  I2N-J2N)  +  ■  •  • 

+  (/ l  -irv- 1*/ 1  -t.v-at  +  .  .  +  /r-j/rll  + 

+  2t(/i//v)  +  1/n+  1/ 2A/ 1  +  .  +  l/i-(N-i|/r)ll 


(B2) 


From  the  definition  of  the  variance  <x2  as  given  in  Eq. 
(2),  the  first  set  of  terms  in  parentheses  on  the  right- 
hand  side  of  Eq.  (B2)  is  just  IV2.  Within  each  curly 
bracket  the  terms  in  parentheses  are  cross-product 
terms  obtained  from  the  corresponding  parentheses  in 
Eq.  (Bl),  and  the  curly  brackets  have  been  arranged 
such  that  the  sets  of  terms  in  the  jth  bracket  are  related 
to  the  jth  correlation  coefficient  pj  [see  Eq.  (3)].  In 
general,  the  relationship  would  be  equal  to  (T  —  j)o2pj, 
except  that  the  fractional  number  j/N  of  the  terms  are 
missing.  These  arise  from  the  absence  of  cross-product 
terms  between  parentheses;  i.e.,  for  j  =  1,  the  terms 
f  n!n+  1,  IznI->n+  1,  etc.  are  missing.  If  we  assume  that 
the  missing  interpolated  terms  have  the  same  average 
value  as  the  terms  present,  Eq.  (B2)  can  be  rewritten 
as 


the  agreement  between  Eq.  (B4)  and  the  experimental 
results  given  in  Tables  I  and  II. 

It  can  be  shown,  assuming  pj  =  constant,  that  the 
fractional  error  introduced  by  ignoring  the  term  (1  — 
j/D  in  Eq.  (B4)  is  approximately  equal  to  jV/3I'.  In  the 
experiments  considered  in  this  paper,  I'  =  12,288; 
therefore,  for  small  N  the  fractional  error  is  negligibly 
small.  For  the  worst  cases  considered,  N  =  1024  and 
2048,  the  fractional  errors  would  be  slightly  <3%  and 
6%,  respectively.  However,  these  represent  an  over¬ 
estimate  of  the  errors,  since  p;  is  not  a  constant  but, 
rather,  a  decreasing  function  of  j.  Therefore,  to  a  good 
approximation,  Eq.  (B4)  can  be  further  simplified  to  the 
form 

,  <t2|  n-i  | 

4»-|l  +  2  L  0-i/N)pyj  •  <B5) 


Appendix  C.  Running  Averaging 

The  variance  of  the  I1  lidar  returns  for  a  running  av¬ 
erage  over  N  returns  is  defined  by 
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+  2 III  1/3  +  lil 4  +  ...  +  Ip-il /  )  +  Uih  +  . . .  +  / r— .■»/ J-— j )  +  .  +  Un-iIn  +  +  / r — <iv - 1 »/ r-iN-aiit  I 
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The  set  of  terms  in  the  first  parentheses  of  the  first 
curly  bracket  on  the  right-hand  side  of  Eq.  (C2)  is  equal 
to  T a2,  and  each  succeeding  set  in  the  bracket  is  similar 
but  with  the  first  and  last  term  of  the  preceding  set 
eliminated.  The  first  set  of  terms  in  the jth  succeeding 
curly  bracket  is  similarly  related  to  the  correlation 
function  pj,  being  equal  to  (T  —  7)  o2pj,  with  succeeding 
sets  again  dropping  the  outer  terms.  As  in  Appendix 


B,  we  assume  that  the  missing  terms  have  the  same 
average  value  as  the  terms  present.  In  Appendix  B  the 
missing  terms  were  interpolated  ones;  here  the  situation 
is  reversed,  with  the  central  terms  present  and  the 
missing  terms  obtained  by  extrapolation  at  both  ends. 
With  the  equal  weighting  assumption.  Eq.  (C2)  be¬ 
comes 
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+  . . .  +  2pN-2<r2 lir  -  (N  -  2)1  +  <r  -  N)l  +  2pN-,ff2 tr  -  (N  - 1)11. 

Since  T  -  2(N  -  1)  +  k  =  (r  -  k)  -  2[N  -  (A  +  D], 
we  obtain 
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Evaluating  the  summations  yields 
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or, 

-  —  [i  +2  E  0  ->/Wp>|  •  (C7) 

The  above  relationship  between  the  correlation 
coefficient  and  the  pulse-averaged  variance  for  a  run¬ 
ning  average  differs  from  the  relationship  obtained  for 
segmental  averaging,  which  is  given  in  Eq.  (B4). 
However,  as  was  discussed  in  Appendix  B,  the  differ¬ 
ence  is  quantitatively  small,  and  the  approximate  form 
for  segmental  averaging  as  given  in  Eq.  (B5)  is  identical 
with  Eq.  (C7). 


Appendix  D.  Analytic  Correlation  Functions 

In  this  Appendix  we  investigate  the  pulse-averaged 
variance  of  the  mean  as  a  function  of  the  number  of 
pulses  averaged  for  three  different  correlation  coeffi¬ 
cients,  pj,  which  are  analytic  functions  of  j.  These 
analytic  functions  are  used  in  Eq.  (5)  to  calculate 
(1 )  Pj  =  constant  =  p: 


N- 1 

P  E  (1  -j/N)‘ 

i-i 


(N-  1) 


— 11  +  PIN 
N 

For  pN  »  1,  a%  *  a2  p. 

(2)  pj  =  exp(->«): 
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